Critical sound attenuation in a diluted Ising system 
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The field-theoretic description of dynamical critical effects of the influence of disor- 
der on acoustic anomalies near the temperature of the second-order phase transition 
is considered for three-dimensional Ising-like systems. Calculations of the sound at- 
tenuation in pure and dilute Ising-like systems near the critical point are presented. 
The dynamical scaling function for the critical attenuation coefficient is calculated. 
The influence of quenched disorder on the asymptotic behaviour of the critical ul- 
trasonic anomalies is discussed. 

PACS numbers: 64.60.Ak, 64.60.Fr, 64.60.Ht, 43.35. Cg 

The progress achieved in the understanding of critical phenomena has largely been due to 
theoretical and experimental works devoted to studying the critical dynamics of condensed 
media. We have seen from experiments (Fig. ^) [lj that for a solid an anomalous peak of 
ultrasonic attenuation is observable in the vicinity of the critical point. The critical anoma- 
lies exhibited by sound attenuation have long been recognized as important in the study of 
dynamical critical phenomena. Ultrasonic methods permit simultaneous measurements of 
both static and dynamic properties. Measurements of sound velocities give information on 
the equilibrium properties, while measurements of the sound attenuation yield information 
on the dynamic properties ofmaterials. The main difficulty in the theoretical discussion of 
the critical propagation of sound waves consists in the estimation of the four-spin correlation 
function. The method which is based on the representation of the four-spin correlation func- 
tion through two-spin correlation functions by means of decoupling leads to overestimated 
values of the critical fluctuations. 

There are a lot of theories and phenomenological descriptions 
ultrasonic anomalies in solids with good agreement with experiments 
materials and crystals have many structural defects and it is worth considering the influence 
of such defects or disorder on the dynamical process of sound propagation in solid media. 




* E-mail : prudnikp@univer.omsk.su 
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FIG. 1: Ultrasonic anomaly for experimental attenuation in FeF2 
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Structural disorder and the presence of impurities or other defects play an important role 
in real materials and physical systems. They may induce new phase transition types of and 
universality classes and modify the d yna mic transport properties. 

According to the Harris criterion j 1 31] - critical behaviour of Ising systems is changed by 
presence of a weak quenched disorder. The problem of the influence of disorder on the critical 
sound propagation in Ising-like systems has been discussed in 14] with use of a e-expansion in 
the lowest order of approximation. However, our pilot analysis of this phenomenon showed 
that in 14| some diagrams which are needed for a correct description of the influence of 
the disorder were not considered. Furthermore, our numerous investigations of pure and 
disordered systems performed in the two-loop and higher orders of the approximation for 
the three-dimensional system directly, together with the use of methods of series summation, 
show that the predictions made in the lowest order of the approximation, especiall y o n the 
basis of the e - expansion, can differ strongly from the real critical behaviour jisi 
Therefore, the results from [14( must be reconsidered with the use of a more accurate field- 
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theoretic approach in the higher orders of the approximation. 

In this paper we have realized the correct field-theoretic description of dynamical effects 
of the influence disorder on acoustic anomalies near the temperature of the second-order 
phase transition for the three-dimensional Ising-like systems in the two-loop approximation. 

In our description, we extended the model of the phase transition in a disordered system 
with a coupling between nonfluctuating variables [17[ to the case (physically important for 
structural phase transitions) of a compressible three-dimensional Ising model with frozen-in 
lattice defects that is considered in using the renormalization-group method in the two-loop 
approximation. 

The interaction of the order parameter with elastic deformations plays a significant role in 
the critical behaviour of the compressible system. It was shown for the first time in [IS] that 
the critical behaviour of a system with elastic degrees of freedom is unstable with respect to 
the connection of the order parameter with acoustic modes and a first-order phase transition 
is realized. However, the conclusions of Q| are only valid at low pressures. It was shown 
in that in the range of high pressures, beginning from a threshold value of pressure, the 
deformational effects induced by the external pressure lead to a change in type of the phase 
transition. 

The Hamiltonian of a disordered Ising model with allowance for elastic degrees of freedom 
may be specified as 

H = H e i + + H int + H imp , (1) 

consisting of four contributions. 

The elastic part is determined from 



H, = \ I d d x 
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where u a/ 3(x) are components of the strain tensor and are the elastic moduli. 
H op is a magnetic part in the appropriate Ginzburg-Landau form: 



/ - < t^t S 2 + ~ (VS 1 ) 2 + ^u S* 



Hrm = / d d X 



(3) 



where S(x) is the Ising field variable which is associated with the spin order parameter, uq is 
a positive constant and r ~ (T — T 0c ) /T 0c with the mean-field phase transition temperature 
Toe- 
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The term H int describes the spin-elastic interaction 



int 



d d x 



go 



(4) 



which is bilinear in the spin order parameter and linear in deformations, go is the bare 
coupling constant. 

The term Hi mp of Hamiltonian determines the influence of disorder and it is considered 



in the following form: 



Himp 



d d x [At(x)S 2 ] + J d d ; 



h(x)J2 



(5) 



where the random Gaussian variables At(x) and h(x) are the local transition temperature 
fluctuations and induced random stress, respectively. Taking the Ar(x) fluctuations into 
account causes additional interaction of order parameter fluctuations over defects and renor- 
malization of the phase transition temperature for disordered systems. Taking into account 
the h(x) fluctuations leads to renormalization of the elastic moduli and the coupling constant 
in the spin-elastic interaction. 

The Fourier transformed variables become 



a/3 



W _|_ v 1/2 u a(3{q) exp (iqx) 



(6) 



with u a p(q) = i/2 [q a up + qpu a ). Then the normal-mode expansion is introduced as u(q) = 
J2x^\(l)Qq,x with the normal coordinate Q q ^\ and polarization vector e\(q). We carry out 
the integration in the partition function with respect to the nondiagonal components of the 
uniform part of the deformation tensor u^l, which are insignificant for the critical behaviour 
of the system in an elastically isotropic medium. 

After all of the transformations ||2fJ the effective Hamiltonian of the system has become 

H=\f d d q (to + q 2 ) S q S-q + lj d d q Ar_ 9 S Ql S q „ qi 

+ ^u J d d q S qi S q2 S q3 S 

-91-122-93 + J d d q qh- q Q q ^\ 
~^wo J d d g {SqS-q) (S q S^ q )-g J d d q qQ- Qi \ S qi S, 
+ a J d d q q 2 Q q ^ Q 



(7) 



J q-qi 



-q,Xi 



where w = 3g 2 /(2V (4C° 12 - C° n )) ; a = (C° n + 4C° 2 - 4C 44 ) /4V 
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The renormalization-group analysis of the critical behaviour of the disordered compressed 
Ising model with the Hamiltonian (JJJ) was carried out in our paper j^J, and different fixed 
points for the Hamiltonian (jZJ and conditions of their stability were determined in the two- 
loop order of the approximation with the use of the Pade-Borel summation technique. It 
was shown that the Ar q fluctuations are relevant for the critical behaviour of the Ising model 
and the replica averaging procedure over Ar q leads to a new vertex of interaction for the 
order parameter fluctuations. It was shown in |20j that random stress connected with h q in 
(J7J) can lead to multicritical behaviour of a system if certain conditions are fulfilled. 

The critical dynamics of the system in the relaxational regime can be described by the 
Langevin equations [2jJ for the spin order parameter S (q) and phonon normal coordinates 

Qx (?) 

dH 

dSi 



S q = -T h 6g + T hi 



-<? 
OH 

Qq,x = q 2 D Q q ,x + V q ,x + h Q , (8) 

where Tq and Dq are the initial kinetic coefficients, £, q (x, t) and rj q (x,t) are Gaussian white 
noises. 

The quantities of interest are the response functions G(q,u) and D(q,cj) of spin and 
deformation variables, respectively. It can be obtained by linearization in correspondent 
fields that 

D{q, u) = 5 [{Q q , w )\ /5h Q = [{QwQ- q ,-Jj\ , (9) 
G(q, lu) = 5 [(S q ,„)} /5h s = [(S q ,„S^)} , (10) 

where (...) denotes averaging over Gaussian white noises, [...] denotes averaging over random 
fields Ar q and h_ q . 

The response functions may be expressed in terms of self-energy parts: 

G-\q,u) = G^(q,uj)+U(q,uj), (11) 
D-\q,uo) = J D - 1 (g, W ) + S(g, W ), 

where the free response functions G (q,u) and D (q,u>) have the forms 

D (q, lo) = 1 /(lj 2 - aq 2 - \ujD q 2 ) , 
G (q,u) = l/(iu/V + (r + q 2 )) ■ 




X(q,u) = 4g* ( > -96g 2 u< Y > + 




+i6 ^C3 > +i6 ^^3 



+ 16g z v 

(a) (b) (c) >-< 



FIG. 2: The diagrammatic representation of the £(</, o>) in the two-loop approximation. Solid line 
(a) corresponds to Go(q,u>), solid line with cross (b) is 2 Go(q, lo)Gq(— q, —uj) /Tq, and vertex with 
dash line (c) corresponds to v = [(At) 2 ] 

The characteristics of the critical sound propagation are defined by means of the response 
function D(q,u). Thus, the coefficient of ultrasonic attenuation is determined through 
imaginary part of H(q,u): 

a(w,T)~wImS(0,w). (12) 



The self-energy part H(q,u>) was obtained by iterative solution |22J of the dynamic equa- 
tion (JHJ) with the effective Hamiltonian ((7j). The diagrammatic representation of £(g, u) in 
the two-loop approximation is presented in Fig. El 

The Feynman diagrams involve momentum integrations on dimension d (in our case 
d = 3). Near the critical point the correlation length £ increases infinitely. When <C A, 
where A is a cutoff in momentum-space integrals (the cutoff A serves to specify the basic 
length scale), the vertex functions are expected to display an asymptotic scaling behaviour 
for wave-numbers g< A. Therefore, one is lead to consider the vertex functions in the limit 
A — > oo. The use of the renormalization-group scheme removes all divergences which arise 
in thermodynamic variables and kinetic coefficients in this limit. 



We have applied the matching method which was introduced for statics in 







23j and gen- 



eralized for critical dynamics in |4|. First, we use the dynamical scaling property of the 
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response function 



D (q, u, t) = etN'D (qe l , (u/T ) e zl , re l/u ) 



(13) 



and then we calculate the right-hand side of this equation for some value 1 = 1*, where the 
arguments do not all vanish simultaneously (2^. The choice of I* is determined by 

2 



It may be argued j^J that condition ()14j) provides an infrared cut-off for all divergent 
values. The particular form of the matching condition (|14j) containing the exponents z and 
v permits an explicit solution for I* 



,2«* 



1. 



(14) 



1 + (y/2f z 



-1/4 



T- u F(y,r). 



(15) 



In (115(1 the abbreviation y = cor zv /T q is introduced and F(y, r) is defined. The value 
of the exponent v = 0.70 was calculated in for the corresponding fixed point. The 
exponent z = 2.1653 was taken from for the 3D disordered Ising model considered for 
purely relaxational model A. The coupling of the order parameter with elastic deformations 
is irrelevant for the dynamics of a disordered system with a negative specific heat exponent 

The response function D (qe l , (lj/Tq) e zl ,re l ^ u ) on the right hand side of ()13() is repre- 
sented by the Dyson equation (jllj) and for the self-energy part we obtain 

Im£(u;) / (a + zu)\ ImSfcue 2 ^ 



CO 



exp / 



v 



toe 



zl 



(16) 



It was shown in later theoretical works that in asymptotic regions the coefficient 

of attenuation can be described using a simple scaling function of the variable y only. The 
experimental investigations performed on the three-dimensional crystals Gd ll| and MnP 



12j confirmed the validity of the concepts of dynamical scaling. 

Thus after renormalization procedure (|16|) we can define a scaling relation in the form 



ImSH /oo =r- a - z y(y), 



(17) 
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FIG. 3: Scaling functions for the critical sound attenuation </>(y) in double logarithmic plotting 
for pure (1) and disordered (2) Ising systems (</>o = 0(0)). 



where <f)(y) is a dynamical scaling function 



2 r n F a/v+1/2u ~ z 



9 *o 



7T 



2 p2z-2/v\ V 2 



1 + 



+ 



1/2' 



+ 



7T" 



*2„,* pa/i 



2p2z-2/v\ 1 / 2 



1/2 



Z-l/l/ 



(4tt) : 



(I/), 



^impvy;j 



(18) 



where g*, u*, and v * are values of vertices in the fixed point of renormalization group 



transformations 



i fimp 



(y) is the numerically calculated contribution of diagrams in £ (a;) 



(Fig. I2J) characterizing the influence of disorder. 

The dynamical scaling function (p(y) is plotted against i/ona double-logarithmic scale 
for pure and disordered systems in Fig. El We have thus seen that the presence of disorder 
is irrelevant for the scaling behaviour in the hydrodynamical region with y <C 1, but it has 
a drastic effect in the critical region with y ^> 1 (T — > T c ). 
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TABLE I: Asymptotical behaviour of the attenuation coefficient in the hydro dynamical and the 
critical regions 



System 


Region 


a(tu, r) 


Pure 


hydrodyn. 


W 2 T -1.37 




critical 


w 1.04 r -0.16 


Disordered hydrodyn. 






critical 


w 1.23 r -0.26 



The anomalous temperature dependence of the calculated attenuation coefficient for pure 
and disordered systems is shown in Fig. UJ Comparison of curves (1) and (2) in Fig. H] clearly 
reveals the strong influence of disorder on the temperature dependence of the attenuation 
coefficient in the vicinity of the critical point. We hope that these theoretical results will 
create a demand for ultrasonic experimental investigations of dilute Ising-like systems, for 
example samples of Fe x Zn 1 _ x F 2 . In producing the figure we used the fact that the presence 
of disorder causes a reduction of the phase transition temperature T c in relation to that of the 
pure system. A model representation of the calculated attenuation coefficient for the pure 
system is shown in Fig.|3]in comparison with experimental data (3) for FeF 2 Qj. Adjustment 
of the experimental data permited us to determine the value of the theoretical parameter 
T and then calculate the attenuation coefficient for the disorder model. We should note 
that the observable differences from the experimental results below T c are explained by the 
contribution of the order parameter relaxation effects to the attenuation, which always occur 
below T c . The relaxation effects are not considered in this paper, but the contribution of 
fluctuations to the attenuation coefficient is relevant over the whole critical range. 

From ()12|) we find that the sound attenuation coefficient obeys the asymptotic scaling 
relation 

a (w, r) ~ uo 2 T~ a ~ vz (j) (y) . (19) 

On the basis of the dependencies ()18j) obtained, we calculated the exponents of asymp- 
totic behaviour of the attenuation coefficient in the hydrodynamical and the critical regions 
(Table We see from the table that in the critical region the anomalies of the attenua- 
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FIG. 4: Thermal dependencies of the ultrasonic attenuation in critical point for pure (1) and 
disordered (2) systems in comparison with experimental results (3) in FeF2 

tion coefficient must be observed in both pure and disordered systems. But for disordered 
systems the anomalies of the temperature and frequency dependences of the attenuation 
coefficient must be stronger than those for pure systems. These conclusions are also shown 
in Fig. |U 

We hope that theoretical estimates of the attenuation coefficient carried out in this paper 
creat a great demand for detection of the anomalies which are found to be induced by 
quenched disorder in ultrasonic experimental investigations of the critical dynamics. 
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